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We study the magnetism of a lattice of coupled tetrahedral spin-1/2 clusters which might be of 
relevance to the tellurate compounds Cu2Te20sX2, with X=C1, Br. Using the flow equation method 
we perform a series expansion in terms of the inter-tetrahedral exchange couplings starting from 
the quadrumer limit. Results will be given for the magnetic instabilities of the quadrumer phase 
and the dispersion of elementary triplet excitations. In limiting cases of our model of one- or two 
dimensional character we show our results to be consistent with findings on previously investigated 
decoupled tetrahedral chains and the Heisenberg model on the 1/5-depleted square lattice. 

PACS numbers: 75.10.Jm, 75.50.Ee, 75.40.-s 



I. INTRODUCTION 

Unconventional magnetism of frustrated spin systems 
has received considerable interest recently. Prominent 
examples are the one-dimeasional (ID) frustrated spin- 
Peierls compound CuGeOatl or the 2D orthogonal spin- 
dimer system SrCu 2 (B03)2 with frustrating inter-dimcr 
couplingsn Apart from dimer-based structures, frus- 
trated spin systems involving triangular or tetrahedral 
units, e.g. the kagome-, checkerboard-, or the pyrochlorc- 
lattices are a focus of current research. In the classical 
limit frustration leads to groiind states with macroscopic 
degeneracy in these systemsau. In the quantum case low- 
lying singlets seem to exist both, on the kagome and the 
checkerboard lattice with no long-range magnetic order in 
the formemEl, and a valencc-bond-crystal (VBC) ground 
state in the latter caseoEl. Analysis of the 3D pyrochlorc 
quantum-magnet remains an open issucO. 

Recently tellurate compounds, Cu2Te205X 2 , with X 
= CI, Br have been found to realize a new class of spin- 
1/2 systems where tetrahedra of Cu 2+ align in tubes 
along the c(z)-direction and are separated by lone-pair 
cations in the a6(ccy)-plandj. Both, the effective dimen- 
sionality of this system as well as the relevant magnetic 
interactions remain a puzzle. Early analysis of thermo- 
dynamic dataallj was based on the OD limit of isolated 
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FIG. 1: a) xy-plane and b) z-axis structure of the 3D tetrahe- 
dral cluster-lattice. Spin-1/2 moments are located on dotted 
vertices. SU(2) type of exchange with strength jo,.. .,6 along 
the links. 



tetrahedral units, i.e. on the exchange pattern of fig. |l| 
with j2,...6 = 0. This resulted in j = 38.5(43)A' and 
ji/jo ~ 1 f° r the Chlorinc(Bromidc) system which has 
been refined recently into jo ~ A7.66K and ji/jo ~ 0.66 
for the Bfpmide system and j 1 /i e -.< 0.66 for the Chlo- 
rine caseO. Raman spectroscopy!^ however, indicates a 
substantial inter-tetrahedral c-axis coupling. This has 
prompted studies of ID tetrahedral spin-chainsli3Ej as 
in fig. Q with ,7*2,3,4,6 = 0. Yet, LDA calculations have 
given evidence of an additional z-axis exchange path 
and transverse inter-chain couplings as shown in a) of 
a magnitude which can not be neglectecO. In fact, spe- 
cific heat data reveals a transition at Tq = 18.2 (11. 4) A 
in the Chlorine(Bromide) system. In the Chlorine case 
the entropy change is consistent with 3D antiferromag- 
netic (AFM) ordering. 

Combining fig. [l] a) and b) a 3D cluster-spin model 
arises, about which very little is known. We believe this 
to be an interesting and highly frustrated magnetic sys- 
tem which deserves to be investigated. Therefore, the 
aim of this work is to shed light onto its excitations and 
possible magnetic instabilities. In addition, our analy- 
sis could be of relevance for the Cu2Te20sX2 system, 
in particular, if additional spectroscopic data becomes 
available. In the remainder of this paper we first discuss 
our method of calculation and then present results on the 
stability and the triplet dynamics. 



II. SERIES EXPANSION 

The Hamiltonian, as read off from fig. [I] can split into 
a bare part Hq and a perturbation Hi 

H = Hq + Hi 

Ho = X>'o(Sii+S 3 i)(S 21 + S4i) 
l 

Hi = [ ji(Sns 3 i + S21S41) 
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where the site of each tetrahedral unit is labeled by 
1 with Sii, i = 1 . . .4 being the spin-1/2 operators cor- 
responding to each tetrahedron and 1 + x(y, z) refers to 
shifts of 1 by one unit cell along the x-r-ij, or z-axis. 

As has been pointed out previouslya the spectrum of 
decoupled tetrahedra, i.e. for J2....6 = 0, is special in as 
such that for ji < jo each tetrahedral ground state is 
a singlet involving all four spins, while at ji > jo it is 
a product of two S'=0-dimers on each of the ji -bonds. 
In turn, at j\ = jo the decoupled local ground states 
are doubly degenerate singlets. This leads to quantum 
criticality in the lattice case_pid_bears the possibility of 
low-lying singlet excitationsQiimij. This may be relevant 
for the Bromide system but in the Chlorine compound j\ 
seems clearly less than jo. Therefore in the remainder of 
this paper we focus on the quadrumer limit of (|l|), defined 
by setting j = 1 and ji,... 6 << 1. 

The spectrum of each quadrumer consists of four 
equidistant levels which can be labeled by spin S and 
a number of local energy quanta q\, c.f. table 1. The un- 
perturbed Hamiltonian Hq which consists of the sum of 
quadrumers displays an equidistant ladder-spectrum la- 
beled by Q = Y^i Qi- The Q = sector is the unperturbed 
ground state |0) of Hq, which is a VBC of quadrumer- 
singlcts. The Q — 1-sector contains local S — 1 single- 
particle excitations of the VBC with qi = 1, where 1 runs 
over the lattice. At Q = 2 the spectrum of Hq has total 
S = 0, 1, or 2 and is of multi-particle nature. For 5 = 
at Q = 2 it comprises of one-particle singlets with q\ = 2 
and two-particle singlets constructed from triplets with 
9i = 9m = 1 and 1 7^ m. The perturbation H\ in (^) 
can be written as a sum of two-site operators T n _k which, 
for each coupling constant jk=i,...e create(destroy) n > 
(n < 0) quanta within the ladder spectrum of Hq. 
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It has been shown reccntlyEHiJli3'E£l that problems of 
type (|^) allow for perturbative analysis using a continu- 
ous unitary transformation generated by the flow equa- 
tion method of WegncrtZL|-The unitarily rotated effective 
Hamiltonian H c g readst-ffia 
00 

H e fi = Hq + Y, <?(m)W roi W„ 



E 

n=1 |m| = n 
Af(m)=0 



■ W„ 



(3) 



where m = (mi . . . m n ) with |m| = n is an n-tuple of 
integers, each in a range of rrii £ {0, ±1, . . . , ±N} and 
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TABLE I: Energy (E, in units of jo), spin (S), and quantum- 
number qi of the quadrumer spectrum. 



W n = J2k=iikT n ,k- In contrast to H of (|l|), H cS con- 
serves the total number of quanta Q. This is evident 
from the constraint M(m) = Ya=i m « = 0- The am- 
plitudes C (m) arc rational numbers computed icpm the 
flow equation methodc3c3. Explicit tabulationEj of the 
T nt k shows that for the Hamiltonian in (|l]) N = 4. 

Q-conservation of H c g leads to a ground state energy 
of Eg = (0\H c h\0). Evaluating this matrix element on 
clusters with periodic boundary conditions, sufficiently 
large not to allow for wrap around at graph-length n 
one can obtain series expansions (SEs) for E g valid to 
0(n) in the thermodynamic limit, i.e. for systems of 
infinite size. Q-conservation also guarantees the Q = 
1-triplcts remain genuine one-particle states. A priori 
single-particle states from sectors with Q > 1 will not 
only disperse via H e g, but can decay into multi-particle 
states. The dispersion of the single-particle excitations is 
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Sim, 00 E° bc are hop- 



where t^i m = {^,lm\H cS \fi, 00) 

ping matrix elements from site (0,0) to site (l,m) for 
a quadrumer excitation /1 inserted into the unperturbed 
ground state. For the thermodynamic limit t^i m has to 
be evaluated on clusters with open boundary conditions 
large enough to embed all linked paths of length n con- 
necting sites (0, 0) to (I, m) at 0{n) of the perturbation. 
E° bc = (0\H cS \0) on the ^,00-cluster. 

Previous applications of this method to spin systems 
were focused on obtaining high-order SEs for one and two 
parameter dimerEEl and quadrumcrQ models in 1 or 2D. In 
the present case, computational constraints related to the 
large number of coupling constants and the 3D nature of 
the model confine the expansion to 4th-order. Moreover, 
explicit display of analytic expressions for the elementary 
triplet dispersion has to be limited to 2nd orderE3. 



III. TRIPLET EXCITATIONS AND MAGNETIC 
INSTABILITIES 



In this section we analyze the triplet dispersion Et (k) 
and the stability of the quadrumer phase against mag- 
netic ordering. We begin by considering the result at 1st 
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FIG. 2: Comparing the triplet dispersion along high- 
symmetry directions of a 2D Brillouin zone as obtained 
from a 6th-op4er plaquette expansion for the 1/5-depleted 
square latticdlil (dashed with error bars) with the 4th-order 
quadrumer expansion for the tetrahedral spin-cluster model 
(solid) at ji = 7, fa = A7, ji — A, and j'3,5,6 = 0, with A = 1 
and 7 = 0.1, 0.3, and 0.5 from top to bottom. (Upper and 
lower edges of error bars refer to 4th- and 5th-order plaquette 
expansion.) 

order in ji e for which we find 

E T (k) = 1 + -(j 4 - 2j 2 )(cos(k x ) + cos(k y j) 

+ y (cos(fc 2 . + k y ) + cos(k x - k y j) 
4 

+ g(j6 - j5)cOs(fc 2 ) . (5) 

Interestingly, this expression depends on three, effective 
exchange coupling-constants only, i.e. a = (j'4 — 2j 2 )/3, 
b = j'3/3, and c = 4(j6 — js)/3. Moreover E^ is indepen- 
dent of ji . In fact we find the dispersion to depend on ji 
starting only at 3rd order. Due to the competition of the 
exchange interactions in ([l]) the effective triplet hopping 
amplitudes a and b in (^]) can be of either sign, even for 
purely AFM j 2 . 4,5,6- 

It is instructive to link (||) to other analytic results 
known from related models. In particular, setting j\ = j 2 
and j3,5,6 = 0, the tetrahedral cluster system of fig. [I] is 
identical to a stack of-JIeisenberg models on the 1/5- 
depleted square latticeO. Bond operator theory (BOT) 
has been applied to this model yielding a triplet disper- 
sion of E T (k) = [1 + 2(j4y 2j 2 )/3(cos(fe c ) + cos(fcy))] 1 / 2 
in the quadrumer phaseEJ. To 1st order this is obvi- 
ously identical to (J5j) with the same setting of parame- 
ters. Similarly, for j'2,3.4,6 = the quadrumer limit of 
fig. [l] maps onto the 'dimerized spin-1 chain sector' of 
the_t£trahedral-chain model of the tellurates studied in 
ref.E3llj. BOT has been applied also to that model, lead- 
ing to E T (k) = [1 - 8j 5 /3cos(£; 2 )] 1/2 . Again, the latter 
is identical to 1st order with (||) with the same choice 
of parameters. While this serves as a consistency check 
for the series expansion we note that BOT, which is ap- 
proximate only, differs from the exact series already at 



2nd order. Additional details on this can be found in the 
appendix. 

To test the quality of our perturbative expansion at 
4th order we compare to the plaquette series-expansion 
of Gel£and and collaborators for the 1 / 5-depleted square 
latticeEil. This is achieved by restricting the parameters 
in (|l]) to those of rcf.Eil, i.e. ji = 7, j 2 = A7, ji = A, 
and j3,5,6 = 0. The plaquette series is a one parameter 
expansion for a 2D model, which allows for expansion 
of £r(k) up to 6th order with respect to A, where the 
unperturbed Hamiltonian incorporates 7 exactly. In fig. 
U we contrast the 4th order results from our six param- 
eter expansion for the 3D tetrahedral spin system with 
the plaquette expansion by considering the triplet dis- 
persion. Despite small deviations which set in upon in- 
creasing 7 the overall agreement is satisfying. Since 7 is 
treated exactly within the plaquette expansion, our 4th- 
order quadrumer series does not coincide with one of the 
edges of the error bars in fig. [| which refer to the 4th- 
and 5th-ordcr plaquette series. 

Next we analyze the stability of the quadrumer phase 
against magnetic ordering by identifying the surface in 
parameter space, closest to ji,...6 = which allows for 
triplet softening, i.e. the occurrence of a wave vector kc 
with £/r(kc) = 0. We emphasize that apart from such 
instabilities, the tetrahedral spin system may exhibit oth- 
ers transitions, as e.g. those related to the lst-order, lo- 
cal quadrumer to dimer-product transition on each of the 
tetrahedra. Here we focus on the triplet softening only. 
To begin, in fig. || we depict the instability surface at 1st 
order as obtained from (^) along with the critical wave 
vectors at which softening occurs. Due to the competing 
interactions several ordering patterns are possible, even 
for AFM couplings only. 

Since the tetrahedral cluster-model contains six ex- 
change coupling parameters we will simplify the stabil- 
ity analysis at 4th order by selecting a subset of them 
only. This selection is based on the effective exchange 
constants at 1st order, i.e. we will focus on the stability 
as a function of a, b, and c setting ji.4,6 = 0. Figure ^ a) 
shows the corresponding instability surface. It has been 




FIG. 3: Left(right) panel: rear(front) view of the stability 
surface of the quadrumer phase at 1st order. Faces are labeled 
by the wave vectors kc of the instabilities and a = (j'4 — 
2j2)/3, b = j 3 /3, and c = 4(j 6 - j 5 )/3. 
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FIG. 4: a) top(bottom) panel: rear(front) view of the stability 
surface of the quadrumer phase at 4th order, b) top(bottom) 
panel: wave vectors of instability for jg < (>)0. Labels refer 
to k c : + =(0,0, tt), □ =(7r,7r,0), A =(ir,0,ir), V =(^,0,0), 
O =(0j0, 0), and X —(n,n,n). Selected points m(p)i i .. .5 refer 
to fig. 

obtained from a numerical search for zeros of the gap of 
the 4th-order triplet-dispersion on a mesh of21x41x2 
points in the a&c-space. For this purpose we have used the 
bare series with no Pade approximations applied. The 
surface is not closed at its extremal extensions in the ab- 
plane, rather the stability analysis has been confined to 
the range of parameters shown in this figure in order to 
comply with the finite range of convergence of the per- 
turbativc result. Only commensurate instability wave- 
vectors have been found within the range of parameters 
investigated. The type of these wave vectors is shown in 
H b). While its shape is deformed with a reduced volume, 
the main features of the 4th-order instability surface are 
still consistent with those at 1st order. We find that the 
additional critical wave-vector types with appear along 
the 'edge-regions', i.e. at 7714,5 and p^ occur within a 
parameter-range of poor convergence of the perturbation 
theory. Therefore, these may be subject to change at 
higher orders. 

Finally we consider the triplet dispersion at critical 
coupling strengths. In fig. |^ we show £r(k) for wave- 
vectors k along high-symmetry directions of the Brillouin 
zone. The exchange parameters have been selected from 
the points Pi,..,5 and mi,. ..5 on the instability surface of 
fig. |i| b). The figure demonstrates a rich variety of k- 
dependencies possible. Since CU2TC2O5CI2 seems to or- 




FIG. 5: Left(right) panel: elementary triplet dispersion in 
the quadrumer phase for j$ > (<)0 at onset of instability. 
Solid(dashed) line refers to 4th(3rd) order series expansion, 
x-axis: path denotes path in Brillouin zone. Subplot labels 
wi(p)i,...5 indicate location of exchange parameters on insta- 
bility surface as in fig. ^ b). Insets refer to exchange param- 
eters at m(p)i,...5. 

der magnetically, inelastic neutron scattering data on the 
tellurates would be interesting in order to choose among 
these dispersion for a set of exchange constants relevant 
to the Chlorine system. To check for the convergence of 
the series expansion fig. || contains both, 3rd- and 4th- 
order results. On those faces of the instability surface 
which appear as continuous deformations of the 1st or- 
der surface of fig. |[ i.e. for p(m) 1^3 the perturbative 
result is well converged. However, within the aforemen- 
tioned edge-regions, i.e. for p(m)4 :5 the convergence is 
insufficient. In particular the critical wave- vectors of the 
instabilities deduced from the 4th-order result within this 
region may be an artifact. This remains to be clarified 
in future analysis. 

To summarize, we have performed a quadrumer series- 
expansion for a three-dimensional tctrahcdral cluster 
spin-system using the flow-equation method. We have 
have shown our results to incorporate and interpolate 
between findings known from previously studied either 
one- or two-dimensional quantum spin-systems which are 
found to be limiting subsets of our model. We have ana- 
lyzed the dispersion of the elementary triplet excitations 
and the stability of the quadrumer phase against mag- 
netic ordering. Future studies will have to contrast this 
type of ordering against other transitions possible in this 
cluster system in order to add more information towards 
a complete quantum phase-diagram. We hope that our 
results may prompt further investigations of the tellurate 
compounds CU2TC2O5X2, in particular inelastic neutron 
scattering studies in order to clarify the relation of the 
cluster spin-model to these materials. 
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In contrast to the lst-order result, a dependence on com- 
bined effective exchange constants only, i.e. (2j2 — j'4) 
and (j5 — j6), is absent. For the case of the 1/5-deplctcd 
square lat tice, i.e. for j\ = j% and J3 = js = je = 0, and 
rewriting ( |AlD in a form which allows for direct compar- 
ison with the BOT of refB3 we get 



APPENDIX A 

To further clarify and connect to other existing ana- 
lytic approaches, in this appendix we also list the result 
to 2nd order in ji t ...e for the triplet dispersion for which 
we find 



J5r(k) = 1 



145jf 31jf 187J2J4 31jf 8jg 
432 864 432 864 27 



x cos(k x ) cos(k y ) + (- y + + -ijfj) ( 2 h - k) 



+ ~72^ ( C0S (M + cos (^)) + — A 3z cos(2k x ) cos(2fc y ) 



1 

5T 



■2 -2 ■ ■ -2 

+ ( - h. - h. + J Ji± + J±_) ( C os(2/U + cos(2fc„)) 
1 27 18 27 108 ;v v ; v v " 

+ -^f33( 2 h - 3j4)(cos(2k x ) cos(k y ) + cos(k x ) cos(2fc. y )) 
+ (-f -fc^ + f)cos(M 

g 

+ 7:J3(j5 - Je) cos(k x ) cos(k y ) cos(k z ) 
4 

- g ( 2 J2 - h){h - k){cos{k x ) + cos(k y )) cos(fc z ) 
-^(j 5 -j 6 ) 2 cos(2fc z ) . (Al) 



E T {k) = 1 



59j 2 2 73j 2 j 4 
144 144 

„-2 



864 



(A2) 



+ (5(j4-2j a )-| + M)(cos(fe) + 

cos(fc y )) 2 



cos(fc y )) 



~18^ 4 ~ 2 ^ 2 ) 2 ( cos (^ 



18 " 

. /4j| 4 j2j4 . 2j|. 2 2 

+ (^--^ + -^)(cos(fc x ) +cos(fc y ) ) . 

this shows the BOT-dispersion, cited after (|5|), to be cor- 
rect only to 1st order. Analogous, for the 'dimerized 
spin.- 1.. chain sector' of the tetrahedral chain stud ied in 
refs.E3tj, i.e. for ,72,3,4,6 = 0, we may rewrite (AL) into 



£r(k) 



, 20i| 
27 

--j 5 cos(fc z 



(A3) 



Againplhc BOT-dispersion is correct to 1st order only. 



In ref.L_j, perturbation theory up to 2nd order has been 
performed using a very different method than pre sen ted 
here. Therefore it is satisfying to realize, that (A3) is 
exactly identical to the corresponding eqn. (9) in section 
II.B.2. ofrefM 
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